UH-511-1075-05 
MADPH-05-1436 



The gluon-fusion uncertainty in Higgs coupling extractions 

Charalampos Anastasiou 

Institute for Theoretical Physics, 
ETH, 8093 Zurich, Switzerland 



Kirill Mclnikov 



in 
o 
o 

(N 



> 

o 
o> 
o 
in 
o 



Department of Physics and Astronomy, University of Hawaii, 
2505 Correa Rd., Honolulu, Hawaii 96822 

Frank Petriello 

Department of Physics, Johns Hopkins University, 
3400 North Charles St., Baltimore, MD 21218 
and 

Department of Physics, 
University of Wisconsin, Madison, WI, 53706 

We point out that the QCD corrections to the gluon-fusion Higgs boson production cross section 
at the LHC are very similar to the corrections to the Higgs decay rate into two giuons. Consequently, 
the ratio of these two quantities has a theoretical uncertainty smaller than the uncertainty in the 
cross section alone by a factor of two. We note that since this ratio is the theoretical input to 
analyses of Higgs coupling extractions at the LHC, the reduced uncertainty should be used; in 
previous studies, the full cross section uncertainty was employed. 
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The benchmark for future new particle searches is the 
search for the Higgs boson at the LHC. It is the only 
particle predicted by the Standard Model (SM) that has 
not been observed, and its discovery will be an important 
milestone for high energy physics. Experiments at LEP 
and SLC established a lower limit on the mass of the 
SM Higgs of ttih > 114 GeV through direct searches pj. 
Indirect constraints from precision electroweak measure- 
ments indicate that a scalar Higgs boson with van < 200 
GeV exists Q • This state will provide a powerful probe of 
the mechanism of electroweak symmetry breaking, which 
has motivated numerous analyses studying which of its 
properties can be measured: whether its CP properties 
can be determined, if its production and decay modes can 
be accurately measured, and a host of other issues |3|]. 

A large effort has been devoted to studying the extrac- 
tion of Higgs boson partial widths at both the LHC and 
a future linear collider |4|. The starting point for these 
analyses is the following relation for each Higgs produc- 
tion and decay channel: 



a p (H) x BR(H -> xx) 



p 



(1) 



where T x denotes the partial width for the decay H — > 
xx, T p is the partial width for the Higgs production 



mode under consideration, and T is the total Higgs width. 
This equation is valid under two assumptions. First, the 
narrow-width approximation is used, which should hold 
up to intermediate Higgs boson masses. Second, it is as- 
sumed that a single, well-defined channel can be assigned 
to the Higgs production process. The observed rate fur- 
nishes a measurement of the product T p T x /T, subject to 
experimental errors and the theoretical uncertainty in the 
normalization factor a p M (H) /Tp M the theoretical input 
and the experimental measurement are clearly separated 
using Eq. ^ The available production modes include 
gluon fusion, weak boson fusion, tiH associated produc- 
tion, and WH, ZH associated production. The accessi- 
ble decay modes in each production channel are: gg — > 

h -» ww,zz,- n ; wbf with h -> ww,zz,^,tt; 

UH production with H -> WW, 77, bb; WH, ZH as- 
sociated production with H — > WW, 77. A combined 
analysis of the many channels available determines the 
various partial widths r^. 

The largest theoretical error entering this analysis is 
the uncertainty in the gluon-fusion production cross sec- 
tion, gg — » H. Next-to-leading order (NLO) QCD cal- 
culations of this cross section revealed corrections reach- 
ing nearly 70% and exhibiting large residual scale depen- 
dences [3 . This stimulated the calculation of the NNLO 
corrections. The inclusive NNLO result for gg — > H was 
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obtained in very recently the leading N 3 LO correc- 
tions to the inclusive cross-section were computed 0. 
The completely differential Higgs production cross sec- 
tion, which allows a realistic study of experimental ac- 
ceptances at NNLO, was calculated in gj. The partial 
width Tgg 1 was also calculated to NNLO 9J. However, 
in spite of this significant effort, the gg — > H channel still 
suffers from large theoretical uncertainties. The residual 
scale dependences in both the cross section and partial 
width are at the 20% level. In the most recent study of 
Higgs couplings at the LHC ^lj > a ±20% theoretical error 
was assigned to the gg — > H channel. The next largest 
error entering the analysis is a 15% uncertainty in the 
tiH production cross section; the remaining theoretical 
errors are under 10%, as are the individual experimental 
systematics. 
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FIG. 1: fiR scale dependence for the Higgs production cross 
section at the LHC, as a function of ran- fiR is varied between 
rrih/2 < [ir < 1m,h, while /if = mu. The LO, NLO, and 
NNLO distributions are shown. 

In this note, we point out that the appropriate uncer- 
tainty in the gg — > H channel which enters the analysis of 
Higgs couplings should instead be ±5%, which is smaller 
by a factor of four. This reduction relies upon the obser- 
vation that the theoretical input for the Higgs coupling 
determination is the ratio o-| 9 M /r|j^. The QCD correc- 
tions to Gg^ 1 and g 1 track each other, and a large por- 
tion of the uncertainty cancels when the ratio is taken. 
We believe that the ±20% error assumed in the analysis 
of i s overly conservative. The NNLO Higgs produc- 
tion itself suffers from an uncertainty of only ±10% [6J, 
and as we will argue later, the imprecise knowledge of 
the gluon distribution function is smaller than the resid- 
ual scale uncertainty. Since this is the largest single er- 
ror entering the analysis of Higgs couplings, its reduction 
may have an important effect on the precisions of the r^. 

We now briefly explain why the theoretical uncertainty 
in the ratio is reduced, and then present a numerical 
proof. The perturbative expansions for the cross section 



and partial width have the forms 

a* g M ~ a 2 MC 2 {l+a s ( f i R )X 1 (^ R ,^ F ) 
+ a 2 s (fi R )X 2 (fiR, Hf) + ■ ■ ■} , 

+ a 2 s ( m )Y 2 ( m ) + ...}. (2) 

fi R and \x f respectively denote the renormalization and 
factorization scales, while C denotes the Wilson coeffi- 
cient obtained from integrating out the top quark. The 
width has no fip dependence, as no initial-state mass- 
factorization is required. An immediate consequence of 
these formulae is that both a^M anc j F^* f are propor- 
tional to the same factor of a 2 (/i R )C 2 . This quantity 
contributes a large renormalization scale dependence to 
both cr?„ M and T^ 1 . The same scale ur can be cho- 
sen in both calculations, which completely removes the 
large renormalization scale uncertainty when the ratio 
cjgg 1 /Fgg 1 is taken. ^Vhe n integrating out the top quark 
to obtain the Wilson coefficient C, the relevant scale is 
the top quark mass m t ; this choice removes all large loga- 
rithms from the expression. This dependence is identical 
in both tjg^ and T?^ , and therefore cancels exactly. 
The observation that /ir is approximately the same for 
both (Tgg 1 and r^'* / in the remaining pieces relies upon 
the fact that the Higgs boson production is dominated 
by partonic threshold. This implies that the partonic 
center-of-mass energy squared s, relevant for the gg — > H 
process, is s ~ m^. Since uih is also the only mass scale 
that enters the decay rate, the kinematic scales for the 
production and decay processes must be similar. This 
indicates that identical {ir should be chosen for both 
processes, and also implies that the {ir dependence of 
the expansion coefficients for the cross section and width 
should be similar, Xi ~ 3^. This further reduces the 
renormalization scale dependence. Since the factoriza- 
tion scale uncertainty turns out to be small, the cancel- 
lation of the renormalization scale dependence leads to a 
reduced theoretical error in the ratio. 

We now present numerical proof of this assertion. We 
estimate the theoretical errors by varying both scales 
within the range mjy/2 < /j,r, /ip < 2m,H- We vary each 
scale in the range vrijijl < Hr,[If < 2m#, and fix the 
other scale to m#. The fiR variation of the cross sec- 
tion is 20% to 25% even at NNLO (see Fig.HJ, while the 
Hp variation is less than 5% at NNLO and is therefore 
negligible compared with hr scale variation. The /ir de- 
pendence is reduced to the 10% level when the ratio is 
taken, as shown in Fig. [21 while the fip dependence is 
unchanged. We note that because of the complete can- 
cellation of the fiR dependence at LO, we do not obtain 
a reliable estimate of the theoreical uncertainty until the 
NNLO corrections are included. We use the program FE- 
HiP |8| to obtain cross section numbers, and we take the 
width results from ||| . We believe that the scale variation 
bands presented here are a good estimate of the theoret- 
ical uncertainty: the NLO and NNLO bands overlap, the 
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fixed order results are in good agreement with those ob- 
tained with threshold resummation |10|. and the leading 
N 3 LO corrections do not significantly change the NNLO 
results 0. 
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FIG. 2: hh scale dependence for the a/T ratio at the LHC, as 
a function of ran- is varied between nih/2 < fiR < 2mu, 
while fiF = mH- The LO, NLO, and NNLO distributions are 
shown. 

In Fig. [31 we estimate the remaining theoretical error in 
O = <yg g M ,^ g M /T s J { by first setting m = fi F = fi, and 
defining AO = \0{ji = 2m H )-0{p = m H /2)|/20 average . 
We can then find the appropriate error band by taking 
±AC In both cases there is a slight cancellation between 
the /ir and dependences; we conservatively round the 
error up and estimate Act w ±10% and A(a/T) w ±5%. 
As claimed, a large portion of the theoretical uncertainty 
cancels in the ratio. 

We also point out that the cancellation of the QCD 
uncertainties between a gg and T gg is expected even after 
mild cuts on the final states are imposed. For example, 
the calculation of |8| shows that the QCD corrections to 
a(gg — > H — > 77) are quite similar to that of the total 
gg — > H cross section even if the standard cuts used 
by ATLAS and CMS for the identification of the two- 
photon signal are applied. We also expect a reduction 
of the theoretical uncertainty in quantities of the type 
a (99 ~* H + jets)/r g9 . However, this reduction should 
be less effective than in the total cross section because 
a(H +jets) depends on additional kinematic factors, such 
as the p± of the jets, that control the magnitude and scale 
dependence of QCD effects. 

Finally, we note that the Higgs boson production cross 
section suffers from imprecise knowledge of parton dis- 
tribution functions (pdfs). The resulting uncertainties 
strongly depend on the production mode and the mass of 
the Higgs boson; they were recently analyzed in 0. For 
gg — > H and a Higgs boson with mass m# ~ 200 GeV, 
the pdf errors lead to an uncertainty in the Higgs bo- 
son production cross section of about ±2.5%, smaller 
than the scale variations studied here and negligible when 



added in quadrature. In contrast to the renormalization 
and factorization scale uncertainties, the pdf error is re- 
ducible, since measurements of Standard Model processes 
at the LHC may lead to more accurate pdf determina- 
tions. For these reasons, we do not consider the pdf un- 
certainty in what follows. 

The ±20% uncertainty on the gluon fusion production 
cross section is the dominant uncertainty in the analysis 
of f° r m H > 150 GeV. While the full re-analysis is 
required to access the impact of the uncertainty reduction 
that we advocate, some simple estimates can be made. 
Consider the determination of T gg from the production 
of TV-bosons in the gluon fusion, gg — ► H — > W + W~ . 
The uncertainty in T gg is given by 




(3) 

where N is the number of events, T — l^ gg ! is the 
theoretical normalization factor and T^y is the decay rate 
for H — > WW. For the purpose of this estimate we 
assume that the errors are not correlated. Following [Tlj |. 
we use SN/N w KT 1 , 8T/T w 1.4 xlO" 1 and ST w /r w 
10" 1 , for m H w 160 GeV. Using 5T/T w 2 x HT 1 , 
we find that the uncertainty on T gg is thirty percent. 
As stressed here, it is more appropriate to use ST/T w 
5 x 10~ 2 ; the uncertainty on T gg then becomes about 
twenty percent. An analysis is underway to determine 
the effect of the error reduction on the extraction of Higgs 
properties at the LHC [T^j. 
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FIG. 3: Theoretical error on a and a /V as a function of Higgs 
mass. The LO, NLO, and NNLO uncertainties are shown. 
The green and red lines at the bottom denote the theoretical 
errors at LO and NLO, respectively. As explained in the 
text, only at NNLO is a reliable error estimate dependence 
obtained. 

We thank T. Plchn and D. Rainwater for many useful 
discussions. The work of K. M. was supported by the US 
Department of Energy under contract DE-FG03-94ER- 
40833 and the Outstanding Junior Investigator Award 



4 



DE-FG03-94ER-40833, and by the Alfred P. Sloan Foun- P420D3620434350, and by the Univeristy of Wisconsin 
dation. The work of F. P. was supported by the National Research Committee with funds granted by the Wiscon- 
Science Foundation under contracts P420D3620414350, sin Alumni Research Foundation. 



[1] R. Barate et al. [ALEPH Collaboration], Phys. Lett. B 

565, 61 (2003) |arXiv:hep-e x/03060 33l. 
[2] [LEP Collaborations] , |arXiv:hep-ex/oil2"uT5] 
[3] For a recent review, see M! Carena and 

H. E. Haber, Prog. P art. Nucl. Phys. 50, 63 (2003) 

|arXiv:hep -ph/0208209 . 
[4] D. Zeppenfeld, R. Kinnuncn, A. Nikitenko 

and E. Richter-Was, Phys. Rev. D 62, 013009 

(2000) |arXiv:h ep-ph/0002036 ; A. Djouadi et al, 
|arXiv:he p-ph/0002258 
[5] S. Dawson, Nucl. Phys. B 359, 283 (1991); A. Djouadi, 

M. Spira and P. M. Zerwas, Phys. Lett. B 264, 440 

(1991). 

[6] R. V. Harlander and W. B. Kilgore, Phys. Rev. Lett. 
88, 201801 (2002) arXiv:hep-ph/0201206 ; C. Anasta- 
siou and K. Melnikov, Nucl. Phys. B 646, 220 (2002) 
[arXiv:hep~p h/0207004 ; V. Ravindran, J. Smith and 
W. L. van Neerven, N ucl. Phys. B 665, 325 (2003) 
|arXiv:hep-ph/0302135|. 

[7] S. Moch and A. Vogt, hep-ph/0508265 

[8] C. Anastasiou, K. Melnikov and F. Petriello, Phys. Rev. 



Lett. 93, 262002 (2004) |arXiv:h ep-ph/0409088 ; 

C. Anastasiou, K. Melnikov and F. Petriello, 

|arXiv:hep -ph/0501130 
[9] K. G. Chetyrkin, B. A. Kniehl and M. Steinhauser, Phys. 

Rev. Lett. 79, 353 (1997) arXiv:hep-ph/9705240 . 
[10] S. Catani, D. de Florian, M. Grazzini and P. Nason, 

JHEP 0307, 028 (2003) arXiv:hep-ph/0306211 . 
[11] M. Duhrssen, S. Heinemeyer, H. Logan, D. Rainwater, 

G. Weiglein and D. Zeppenfeld, Phys. Rev. D 70, 113009 

(2004) arXiv:hep-ph/0406323 . 
[12] ATLAS collaboration, "Atlas detector and physics 

performance; technical design report.", report 

CERN/LHCC 99-15, ATLAS-TDR-15; 

CMS collaboration, "CMS: The electromagnetic 

calorimeter, technical design report.", report 

CERN/LHCC 97-33, CMS-TDR-4. 
[13] A. Djouadi and S. Ferra g, Phys. Lett. B 586, 345 (2004) 

arXiv:hep-ph/0310209 . 
[14] D. Rainwater, private communication. 



